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We consider cumulant moments (cumulants) of the thrust distribution using predictions of the full 
spectrum for thrust including 0{al) fixed order results, resummation of singular N'^LL logarithmic 
contributions, and a class of leading power corrections in a renormalon-free scheme. From a global 
fit to the first thrust moment we extract the strong coupling and the leading power correction matrix 
element Qi. We obtain as{mz) = 0.1140 ± (0.0004)exp ± (0.0013)hadr ± (0.0007)port, where the 1-a 
uncertainties are experimental, from hadronization (related to fli) and perturbative, respectively, 
and Cli = 0.377 ± (0.044)cxp ± (0.039)pcrt GeV. The n-th thrust cumulants forn > 2 are completely 
insensitive to ili, and therefore a good instrument for extracting information on higher order power 
corrections, n'„/Q", from moment data. We find (^2)^/^ = 0.74 ± (O.ll)exp ± (0.09)pert GeV. 



I. INTRODUCTION 

The process e+e^ — > jets plays an important role in 
precise determinations of as{mz), as well as for prob- 
ing the nonperturbativc dynamics of hadronization in jet 
production. A wealth of high precision data with percent 
level uncertainties, is available for jet production in e+e~ 
collisions at the Z-pole, Q = mz, and with somewhat 
larger uncertainties at both lower and higher energies Q. 
For a review of classic work on as(™z) determinations 
using event shapes and other jet observables, the reader is 
referred to [T] . Accurate predictions for event shapes are 
now available which include ©(a^) corrections a 
next-to-next-to-next-to-leading-log (N^LL) resummation 
of large logarithms [SI [7j , and a high precision method 
developed for simultaneously incorporating field theory 
matrix elements for the power corrections [8, . 

The majority of fits for as{rnz) from event shapes e 
make use of cross section distributions dcr/de, in a region 
where nonperturbativc effects enter as power corrections 
in \/Q and the theoretical description is the most ac- 
curate. In our recent analysis [5] for the event-shape 
variable thrust t — 1 — T , 

T ^ max^ ' , (1 

\Pi\ 

we obtained a precise determination of as{mz). Our 
theoretical description is based on Soft-CoUinear Effec- 
tive Theory (SCET) [rUHTi] . and has several advanced 
features, such as: 

1. Matrix elements and nonsingular terms at order 
using results from [2]. Non-logarithmic terms in 
the hard function are included at order as well. 

2. Resummation of the singular logarithmic terms to 
all orders in up to N'^LL order. 

3. Profile functions (r-dependent scales /is, -R, 
fins) that correctly treat the peak region and ac- 



count for the multijet boundary condition to ensure 
that predictions converge properly into the known 
fixed order result in the multijet endpoint region. 
They allow an accurate theoretical description over 
the entire range r S [0, 0.5]. 

4. Description of nonperturbativc effects with field 
theory and a fit to a single nonperturbativc matrix 
element of Wilson lines in the tail region (where 
power corrections are described by an OPE). 

5. Definition of f2i in a more stable Rgap scheme 
^T5[ [TC] rather than in MS. This ensures fii and 
the perturbative cross section are free of ©(Aqcd) 
renormalon ambiguities. An RGE is used to sum 
large logarithms in the perturbative renormalon 
subtractions [T71|TH]. The fit gives fii with an ac- 
curacy of 16%. 

6. QED final state corrections at 0{a) and NNLL 
(counting a ^ a^); bottom mass corrections are 
included using a factorization theorem with log re- 
summation; 0{a'^) axial-singlct terms arising from 
the large top-bottom mass splitting are included as 
well. 

A two-parameter global fit in the tail of the thrust dis- 
tribution gives [8] asimz) = 0.1135 ± (0.0002)oxp ± 
(0.0005)hadr ± (0.0009)po,.t as well as f^i = 0.323 ± 
(0.009)exp ± (0.013)n,± (0.020),^(„,,) ± (0.045)pert GeV 
where Qi = fii(i?A,MA) is defined in the Rgap scheme 
at the scales Ra — /^a — 2 GeV. For as the three uncer- 
tainties are the experimental uncertainty, hadronization 
uncertainty coming mainly from the determination of fii, 
and the perturbative theoretical uncertainty. This result 
for as is one of the most precise in the literature. It is 
also one of the lowest, being 3.9 cr away from the 2009 
world average [12] and 4.0 cr from the 2011 world aver- 
age [20]. For a detailed discussion of as{mz) determina- 
tions see Ref. [21]. The small value of as{mz) is directly 
connected to the non-negligible correction from fii [H], 
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whose fit value is of natural size fii ~ Aqcd- Given the 
discrepancy, further tests of the theoretical predictions 
for event shapes are warranted. In this paper we will 
do so using experimental moments involving the thrust 
variable. 

The property of the N^LL + C'(ag) predictions for 
dcr/dr in Ref. [5] that we will exploit is that they are 
valid in both the dijet and tail regions, where singular and 
large logarithmic terms in need of resummation arise, and 
in the multijet region, where fixed order results without 
log resummation should be used. That is, they are valid 
for all values of r (an improvement over earlier results 
at this order). Important ingredients are: the inclusion 
of the nonsingular terms, important away from the peak 
region; the use of profile functions that turn off resum- 
mation in the far-tail region; and the inclusion of a soft 
function, which is necessary to describe the peak in the 
dijet region, where nonperturbative effects are 0(1). 

We will use the full t range results to analyze moments 
Mn of the thrust distribution in e+e^ — > jets, 

M„ - - / dr r" ^ . (2) 

Unlike for tail fits, the entire physical r range contributes, 
providing sensitivity to a different region of the spec- 
trum. Experimental results are available for many val- 
ues of Q, and the analysis of systematic uncertainties is 
to a large extent independent from that for the binned 
distributions. Thus the outcome for a fit of data for the 
first moment Mi to Usirnz) and fii serves as an im- 
portant cross check of the results obtained in Ref. [8]. 
The Mn moments are also not sensitive to large loga- 
rithms, and hence provide a non-trivial check on whether 
the N'^LL + 0(0;^) full spectrum results, which contain a 
summation of logarithms of r with a substantial numeri- 
cal effect for small r values, can reproduce this property. 
We explore this issue both for central values and for the- 
ory uncertainty estimates. 

The second purpose of this work is to discuss the struc- 
ture of higher order power corrections in thrust moments. 
We find that cumulant moments M'^ (cumulants) are very 
useful, since they allow for a cleaner separation of the 
subleading nonperturbative matrix elements compared to 
the Mn moments of Eq. (pi) . Cumulants include the vari- 
ance M'2 and skewness M^, and we will consider the first 
five: 

M[ = Ml , (3) 
M'^ = M2 - Ml , 
M3 = M3 - 3 M 2 Ml -I- 2 Ml , 
M'^ = M4 - 4 M3 Ml - 3 M| -I- 12 M2 Ml -6M^ , 
= M5 - 5 M4 Ml - 10 Ma M2 + 20 M3 Ml 
+ 30 M|Mi - 60 Ml M2 + 24 Mf . 

In the leading order thrust factorization theorem the 
power correction matrix elements for the moments Af„ 



are called 17^ while for the cumulants M^ they are called 
il'n^. (The n'n^ are also related to the by Eq. ^ 
with Mn ^ ^n ■ ) In particular, the invariance of the 
cumulants to shifts in t implies that the M^>2 moments 
are completely insensitive to the leading thrust power 
correction parameter fli, and hence can provide non- 
trivial information on the higher order power corrections 
which enter as fl'n/Q" and as power corrections 

from terms beyond the leading factorization theorem. In 
contrast, for each M„>2 there is a term ^ as^i/Q that 
for larger Qs dominates over the Vlm/Q"^ terms. "'^ 



A. Review of Experiments and Earlier Literature 

Dedicated experimental analyses of thrust moments 
have been reported by various experiments: JADE [22] 
measured the first moment at Q = 35, 44 GeV, and 
in [25] reported measurements of the first five moments 
at Q = 14, 22, 34.6, 35, 38.3, 43.8 GeV; OPAL [24] mea- 
sured the first five moments at Q = 91, 133, 177, 197 
GeV, and there is an additional measurement of the first 
moment at Q = 161 GeV [ISj; ALEPH [H] measured the 
first moment at Q = 91.2, 133, 161, 172, 183, 189, 196, 
200, 206 GcV; DELPHI [27] has measurements of the 
first moment at Q = 45.2, 66, 76.3 GeV, measurements 
of the first three moments at Q = 183, 189, 192, 196, 200, 
202, 205, 207 GeV 28], and at Q = 91.2, 133, 161, 172, 
183 GeV [55]; L3 [3D] measured the first two moments at 
Q = 91.2 GeV and other center of mass energies which 
are superseded by the ones in [3T] at Q = 41.4, 55.3, 65.4, 
75.7, 82.3, 85.1, 130.1, 136.1, 161.3, 172.3, 182.8, 188.6, 
194.4, 200.2, 206.2 GeV; TASSO measured the first mo- 
ment at Q = 14, 22, 35, 44 GeV [21]; and AMY mea- 
sured the first moment at Q = 55.2 GeV [S^. Finally, 
the variance and skewness have been explicitly measured 
by DELPHI [H] at Q = 133, 161, 172, 183 GeV; and 
OPAL [551 at Q = 161 GeV. All of the experimental mo- 
ments will be used in our fits, with the exception of the 
results in Ref. [53] and data with Q < 22 GeV where our 
treatment of &-quark mass effects may not suffice. 

In principle the JADE results in Ref. [53] supersede 
the earlier analysis of this data reported in Ref. [55] . In 
the more recent analysis the contribution of primary hh 
events has been subtracted using Monte Carlo genera- 
tors.^ Since the theoretical precision of these generators 
is significantly worse than our N'^LL -I- ©(af) treatment 
of massless quark effects and our NNLL + 0(q;s) treat- 
ment of TTifc-dependent corrections, it is not clear how 
our code should be modified consistently to account for 



^ The cumulants begin to differ for n > 4 from the so-called central 
moments, ({t — Mi)"). Both cumulants and central moments 
are shift independent, but the cumulants are slightly preferred 
because they are only sensitive to a single moment of the leading 
order soft function in the thrust factorization theorem. 

^ We thank C. Pahl for clarifying precisely how this was done. 
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these subtractions. Comparing the old versus new JADE 
data at Q = 44GeV one finds Mi = 0.0860 ± 0.0014 
versus Mi = 0.0807 ± 0.0016. This corresponds to a 
3.4 CT change assuming 100% correlated uncertainties (or 
a 2.6 CT change with uncorrelated uncertainties). In our 
analysis we find that the older JADE data provides more 
consistent results when employed in a combined fit with 
data from the other experiments (related to smaller 
values). For this reason our default dataset incorporates 
only the older JADE moment data. We will report on the 
change that would be induced by using the new JADE 
data if we simply ignore the fact that the hh events were 
removed. 

Event shape moments have also been extensively stud- 
ied in the theoretical literature. The C(aj) QCD cor- 
rections for event shape moments have been calculated 
in Ref. jSH 135]. The leading A/Q power correction to 
the first moment of event shape distributions were first 
studied in [36H39] often with the study of renormalons 
(see [3D], and [H] for a review). Ref. [H] made a renor- 
malon analysis of the second moment of the thrust distri- 
bution, finding that the leading renormalon contribution 
is not \/Q^ but rather \/Q^. Hadronization effects have 
also been frequently considered in the framework of the 
dispersive model for the strong coupling US] HI] ^ . 
In this approach an IR cutoff is introduced and the 
strong coupling constant below the scale /i/ is replaced by 
an effective coupling aoff such that perturbative infrared 
effects coming from scales below /x/ are subtracted. In 
the dispersive model the term /i/ao is the analog of the 
QCD matrix element Sli that is derived from the oper- 
ator product expansion (OPE). Since in the dispersive 
model there is only one nonperturbative parameter, it 
does not contain analogs of the independent nonpertur- 
bative QCD matrix elements ^n>2 of the operator prod- 
uct expansion. Thus measurements of f^^>2 can be used 
as a test for additional nonperturbative physics that go 
beyond this framework. 

The dispersive model has been used in Refs. 
together with 0{a^g) fixed order results to analyze event 
shape moments, fitting simultaneously to as{mz) and 
OfQ. Recently these analyses have been extended to C(a;^) 
in Ref. [31] , based on code for n/ = 5 massless quark fla- 
vors, using data from [231 El] and fitting to the first five 
moments for several event-shape variables. Our numer- 
ical analysis only considers thrust moments, but with a 
global dataset from all available experiments. A detailed 
comparison with Ref. ^48, will be made at appropriate 
points in the paper. Theoretically our analysis goes be- 
yond their work by using a formalism that has no large 
logarithms in the renormalon subtraction, includes the 
analog of the "Milan factor" [331 SS] in our framework 



^ Another approach to hadronization corrections to moments of 
event shapes distributions based on renormalons is that of Gardi 
and Grunberg 1451. 



at 0{a^) (one higher order than [IS]), and incorporates 
higher order power corrections beyond the leading shift 
from r^i. We also test the effect of including resumma- 
tion. 



B. Outline 

This article is organized as follows: We start out by 
defining moments and cumulants of distributions, and 
their respective generating functions in Sec. [llj where we 
also discuss the leading and subleading power corrections 
of thrust moments in an OPE framework. In Sec. IIIII we 
present and discuss our main results for as{mz) from 



fits to the first thrust moment Mi. In Sec. VI we analyze 
higher moments M„>2. Sec. VII contains an analysis 



of subleading power corrections from fits to cumulants 
M^>2 obtained from the moment data. Our conclusions 



are presented in Sec. VIII 



II. FORMALISM 

A. Various Moments of a Distribution 

The moments of a probability distribution function 
p{k) are given by 



M„ = (fc") = / dfcp(fc)fc". 



(4) 



The characteristic function is the generator of these mo- 
ments and is defined as the Fourier transform 



Piy) 



-iky\ 



dkp{k)e-''^y ^^M„, (5) 



with Mq = 1. The logarithm of p{y) generates the cumu- 
lants (or connected moments) M^ of the distribution 



71=1 



(6) 



and is called the cumulant generating function. For n > 2 
the cumulants have the property of being invariant under 
shifts of the distribution. Replacing p{k) — >■ p(k — fco) 
takes p{y) e~'y'''> p{y), which shifts A/( M[ + fco 
while leaving all M^>2 unchanged. Writing 



22 i^r! = exp 



N=0 



m 



E 

3 = 1 



Rl 



(7) 



i=l R=0 

one can derive an all-n relation between moments and 
cumulants of a distribution: 

M. = mY.Y[^^. (8) 



^ Ki,! 
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Here the Kij are non-negative integers which determine a 
partition of the integer N through J '^ij ~ ■• ^^'^ 

p{N) is the the number of unique partitions of A^. (A 
partition of TV is a set of integers which sum to N. Here 
the number of times the value j appears as a part 
in the I'th partition, and corresponds to R in Eq. ([t]). ) 
As an example we quote the relation for = 4 which 
has five partitions, p(4) = 5, giving 



Ma 



--M^ + 4 M[ + 3 Aff + 6 + M[* 



(9) 



In the fourth partition, i = 4, we have kh — 2, K42 — 1, 
and K43 = K44 = 0, and the factorials give the prefactor 
of 6. Eq. ([s]) gives the moments Mi in terms of the cumu- 
lants Ml, and these relations can be inverted to yield the 
formulas quoted for the cumulants in Eq. ([s]). > is 
the well known variance of the distribution. Higher order 
cumulants can be positive or negative. The skewness of 
the distribution M3 provides a measure of its asymmetry, 
and we expect Mg > for thrust with its long tail to the 
right of the peak. The kurtosis M4 provides a measure of 
the "peakedness" of the distribution, where M4 > for 
a sharper peak than a Gaussian.'^ 

The shift independence of the cumulants Af^ make 
them an ideal basis for studying event shape moments. 
In particular, since the leading ©(Aqcd/Q) power cor- 
rection acts similar to a shift to the event shape distri- 
bution fMl fi5l [5DH5^ . we can anticipate that M^>2 will 
be more sensitive to higher order power corrections. We 
will quantify this statement in the next section by using 
factorization for the thrust distribution to derive factor- 
ization formulae for the thrust cumulants in the form of 
an operator product expansion. 



B. Thrust moments 



We will first make use of the leading order factor- 
ization theorem, da/dr = / dp{da/dT){T — p/Q)Fr{p), 
which is valid for all t. It separates perturbative da/dr 
and nonperturbative Ft-{p) contributions to all orders in 
Us and Aqcd/('9t), but is only valid at leading order 
in Aqcd/Q- For this factorization theorem we follow 
Ref. [5| (except that here we denote the nonperturbative 
soft function by Fr).^ We will then extend our analysis 
to parameterize corrections to all orders in Aqcd/Q- 



Taking moments of the leading order dcr/dr gives^ 
'•'3^ 1 da 



d.."/ ^P-,£{r-^)FAp) (10) 



^ l^drdpe[r„,-r-^ 



p \ " 1 do- , , 



E 



';)(^) 



where a is the perturbative total hadronic cross section 
and all hatted quantities are perturbative. In the last 



line of Eq. ( 10 ) we used 6{Tm — T — p/Q) = 9{Tra — r) [1 — 

^{P/Q - Tm)^ 0{T,n - p/Q) 0{p/Q + T - T,^)] tO obtaiu 

the three terms. In Eq. ( [lO| the term in square brackets 
is our desired result containing the perturbative Af„ and 
nonperturbative r2„ moments 



M„ = 



dr r" --^{t), 
cr dr 



Afo = 1 , (11) 

a. =^ dp (|)'V,(p), n^^i. 



The small "error" terms in Eq. ( 10 1 are given by 



4-^|:C;)(^)"-v,£*(|)"-v,„), 

For the contribution E^'^ the p-integral is smaller than 
2Q-30 £qj. g^j^y Q -[-Qj. -j-j^g jjj,g^ moments, and hence 

Ei^^ ~ 0. This occurs because Fr{p) falls off exponen- 
tially for p > 2 ill ^ 2Aqcd milS^, and hence values 
p> QTm = Q/'i, are already far out on the exponential 
tail. The E), ' term gives a small contribution because 
the integral is suppressed by either Fj- or dcr/dr: near 
the endpoint r ~ — 2 Aqcd IQ the p-integration is not 
restricted and Fr{p) ^ 1, but da/dr is highly suppressed. 
For smaller r the p-integration is restricted and the expo- 
nential tail of Ft-(p) suppresses the contribution. We have 
checked numerically that at Q = 91.2 GeV [Q = 35 GeV], 
for the first moment the relative contribution of E\ 



compared to the term in square brackets in Eq. ( 10 1 is 



0(10"^) [O(IO-S)], while for the fifth moment E^^ it 
is 0(10^®) [0(10^"*) ]. This suppression does not rely on 
the model used for Ft{p). Thus En can also be safely 
neglected. 



* The cumulants of a Gaussian are all zero for n > 2, and the 
cumulants of a delta function are all zero for n > 1. 
Earlier discussions of shape functions for thrust can be found in 
Refs. I53ll54) . 



® This manipulation is valid when the renormalization scales of 
the jet and soft function which implement resummation are in = 
~ P/Q)y rather than the more standard l^iir) used in [8]. 
Both choices are perturbatively valid, and we have checked that 
the difference is 0.4% for Mi, rising to 0.8% for M5, and hence 
is always well within the perturbative uncertainty. 



5 



Within the theoretical precision we conclude that the 
leading factorization theorem for the distribution yields 
an operator product expansion that separates perturba- 
tive and nonperturbative corrections in the moments 



M„ = 



Ml ^In-e ■ 



(13) 



For M„ the terms that numerically dominate are M„ 
and Mn-i^i/Q. However for the cumulants there 



are cancellations, and Eq. (13) does not suffice due to 
our neglect so far of (Aqcd /QT"' suppressed terms in the 
factorization expression for the thrust distribution. 

To rectify this we parameterize the (Aqcd/Q)"' power 
corrections by a series of power suppressed nonpertur- 
bative soft functions, A^~^Ftj{p/A) 



^^QCD- 



Here 



A ^i^T-o(p/A) = Ft{p) is the leading soft function from 
Eq. (10). We introduced the parameter A = 400 MeV ^ 



Aqcd to track the dimension of these subleading soft 
functions. This parameterization is motivated by the fact 
that subleading factorization results can in principle be 
derived with SCET [56] . and at each order in the power 
expansion will yield new soft function matrix elements. 

Both the factorization analysis and calculation of cu- 
mulants is simpler in Fourier space, so we let 



(14) 



^..(^A)./fe--F.,(|), 

and likewise for the leading power partonic cross section 
da/dTlr) — > (To(2/)- The factorization-based formula for 
thrust is then 



(15) 



where c^jyoiy) accounts for perturbative corrections in 
the (Aqcd/Q)"' power correction. The j = term 
is equivalent to the result used in Eq. (10), Fr{p) — 



AFrfl{p/ A), and the normalization condition for the lead- 
ing nonperturbative soft function is Ft q{z = 0) = 1. The 
terms in Eq. ( 15 ) beyond j = are schematic since in 



reality they may involve convolutions in more variables 
in the nonperturbative soft functions (as observed in the 
subleading 6 — )■ s 7 factorization theorem results [SSH5T]). 
Nevertheless the scaling is correct, and Eq. (15) will suf- 



fice for our analysis where we only seek to classify how 
various power corrections could enter higher moments or 
cumulants. 

The identities a{y = 0)/<t = 1 and <To{y = 0)/o- = 1 
together with Eq. (flSj) imply 



F,,,(2/ = 0) = 0, forj>l. (16) 
Using the Fourier-space cross section the moments are 



1 



dy" 



o-(2/) 



(17) 



y=o 



d" 

dy'^ 



^ 00 



00 ^ 



J=0 
n 

E 

i=0 



^Av) ( 



o) ' 



y=o 



which extends the OPE in Eq. ( 13 ) to parameterize the 



(Aqcd / Qy power corrections. Here the perturbative and 
nonperturbative moments are defined as 



d" 

dy" 

■n 

2" dz~^ 



1 



o'i(y) 



A^ Fr,j{zA) 



(18) 



where M„j- is a dimensionless series in as{p) and r2„.j 
Aqcd- order for M„j to exist it is crucial that our 
cTj{y) and its derivatives do not contain ln(j/) dependence 
in the y — limit at any order in as- In r-space the 
perturbative coefficients have support over a finite range, 
r e [0,1/2], and 



1/2 



dr e' 



try 



a,{r). 



(19) 



1/2 

Therefore the existence of /q dr (r) implies a well 
defined Taylor series in y under the integrand in Eq. ( 19 ), 



and hence the existence of Mnj. This integral is the 



total perturbative cross section for j = 0. From Eq. ( 16 ) 
we have rioj>o = 0, and furthermore ft„fi 
Mn,0 = M„. 



r2„ and 



For the first moment, Eq. (17) yields 



Ml = Ml 



3=0 



(20) 



where the first two terms are determined by the leading 
order factorization theorem, while the last term identi- 
fies the scaHng of contributions from (Aqcd / Q)^^-' power 



corrections. Two properties of Eq. ( |20[ ) will be relevant 
for our analysis: first, there is no perturbative Wilson 
coefficient for the leading 2ili/Q power correction; and 
second, terms from beyond the leading factorization the- 
orem only enter at 0{Aqqy)/Q^) and beyond. For higher 
order moments, n > 2, we have 



M„ = Mn + AIn-l + ^ ' 



2 n fil l 



(21) 



Next we derive an analogous expression for the n-th 
order cumulants for n > 2, which are generated from 
Fourier space by 



d" 
djT 



In 



(22) 



6 



Eq. (15 1 can be conveniently written as the product of 



1.10 



three terms 
1 



1 



1 



yA 

Q 
Q 



(23) 



yA 



where bars indicate the ratios 

'^o(y)' 



,iy) 



Frfli^y 



(24) 



From Eq. (16 1 we have Frj{x = 0) = for all j > 1. 
Taking the logarithm of Eq. (23) expresses the thrust 



cumulants by the sum of three terms 



ML = ML + 



r oo 



2 

Q 

(y) 



n' 



d" 

dy" ^ 

" k=l 



(-1) 



fc+i 



q) 



yA 



1 k 



(25) 



The first two terms involve the perturbative cumulants 
and the cumulants of the leading nonperturbative 
soft functions fii, 



M' = 



n' = 



d" 

d" 
dz'' 



In-CTo(y) 

(T 



(26) 



z=0 



The third term in Eq. ( 25 ) represents contributions from 



power-suppressed terms that are not contained in the 
leading thrust factorization theorem. These terms start 
at ©(Aqqp/Q^). At this order only Fr^i has to be con- 
sidered. The terms i^r,i>2 do not contribute due to ex- 
plicit powers of Aqcd/Q- Concerning F^ 2j it must be 
hit by at least one derivative because -Ft,2(0) = 0, and 
hence does not contribute as well. Performing the 7i-th 
derivative at y — and keeping only the dominant term 
from the power corrections gives the OPE 



ML = ML + 



2^ 



- nM„ 



2ni, 



Here flu is defined in Eq. (18) 
ficient is 



(27) 

The perturbative coef- 



M,,i = 



, d^ 
dy^ 



^liy) (28) 

J y=0 

and so far unknown. For n = 2 the absence of a 1/Q 
power correction in Eq. ([27]) was discussed in Ref. [S3]. 

The majority of our analysis will focus on Mi where 
terms beyond the leading order factorization theorem are 
power suppressed. For our analysis of M„>2 we consider 
the impact of both as^li/Q corrections, and power cor- 
rections suppressed by more powers of 1/Q. When we an- 
alyze M^>2 will consider both l/Q" and power 
corrections in the fits. 



cr 



1.00- 



0.90 



0.80- 
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A Resummation with Rgap 



1 2 3 
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FIG. 1: Theoretical computations at various orders in per- 
turbation theory for the total hadronic cross section at the 
Z-pole normalized to the Born-level cross section (Tq. Here the 
small blue points correspond to fixed order perturbation the- 
ory, green squares to resummation without renormalon sub- 
tractions, and red triangles to resummation with renormalon 
subtractions. 



III. RESULTS FOR Mi 

In this section we present the main results of our anal- 
ysis, the fits to the first moment of the thrust distribution 
and the determination of as{mz) and Vli. Prior to pre- 
senting our final numbers in Sec. |IIID| we discuss various 
aspects important for their interpretation 



In Sec. ITiri" 

we discuss the role of the log-resummation contained in 
our fit code, the perturbative convergence for different 
kinds of expansion methods, and we illustrate the numer- 
ical impact of power corrections and the renormalon sub- 
traction. We also briefly discuss the degeneracy between 
as{mz) and fii that motivates carrying out global fits 
to data covering a large range of Q values. In Sec. |IIIB| 
we present the outcome of the theory parameter scans, 
on which the estimate of theory uncertainties in our fits 
are based, and show the final results. We also display 
results for the fits at various levels of accuracy. Sec. |I11 C| 
briefly discusses the effects of QED and bottom mass 
corrections. Sec. IIVI shows the results of a fit in which 
renormalon subtractions and power corrections are in- 
cluded, but resummation of logs in the thrust distribu- 
tion is turned off. 

For our moment analysis we use the thrust distribution 
code developed in Ref. [H], where a detailed description 
of the various ingredients may be found. We are able to 
perform fits with different level of accuracy: fixed order 
at 0(aj), resummation of large logarithms to N'^LL accu- 
racy^, power corrections, and subtraction of the leading 
renormalon ambiguity. Recently the complete calcula- 
tion of the 0{a^g) hemisphere soft function has become 



^ Throughout this publication N"LL corresponds to the same or- 
der counting as N"LL' in Ref. [8]. 
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FIG. 2: Theoretical prediction for the first three moments at 
the Z-pole at various orders in perturbation theory. The blue 
circles correspond to fixed order perturbation theory (normal- 
ized with the total hadronic cross section) at 0{as), 0{ag) 
and 0{al), green squares correspond to resummed predictions 
at NLL, NNLL, and N^LL normalized with the total hadronic 
cross section, and red triangles correspond to resummation 
normalized with the norm of the resummed distribution. For 
these plots we use as{mz) ~ 0.ff4. 



available |62H64] , so the code is updated to use the fixed 
parameter S2 = —40.6804 from Refs. [SHIM]- A feature 
of our code is its ability to describe the thrust distribution 
in the whole range of thrust values. This is achieved with 
the introduction of what we call profile functions, which 
are r-dependent factorization scales. In the e+ anni- 
hilation process there are three relevant scales: hard, jet 



and soft, associated to the center of mass energy, the jet 
mass and the energy of soft radiation, respectively. The 
purpose of r-dependent profile functions for these scales 
is to smoothly interpolate between the peak region where 
we must ensure that /i^ > Aqcd, the dijet region where 
the summation of large logs is crucial, and the multijet 
region where regular perturbation theory is appropriate 
to describe the partonic contribution [8j. The major part 
of the higher order perturbative uncertainties are directly 
related to the arbitrariness of the profile functions, and 
are estimated by scanning the space of parameters that 
specify them. For details on the profile functions and 
the parameter scans we refer the reader to App. |Aj We 
note that our distribution code was designed for Q values 
above 22 GeV. 



A. Ingredients 

The theoretical fixed order expression for the thrust 
moments contain no large logarithms, so we might not 
expect that the resummation of logarithms in the thrust 
spectrum will play a role in the numerical analysis. We 
will show that there is nevertheless some benefit in ac- 
counting for the resummation of thrust logarithms. This 
is studied in Figs. [T] and [2j where for Q — mz we compare 
the theoretical value of moments of the thrust distribu- 
tion obtained in fixed order with those obtained includ- 
ing resummation. (The error bars for the fixed order 
expansion arise from varying the renormalization scale /i 
between Q/2 and 2 Q and those for the resummed results 
arise from our theory parameter scan method.) 

In Fig. [1] we show the total hadronic cross section a 
from the fixed order expansion (blue points with small 
uncertainties sitting on the horizontal line) and deter- 
mined from the integral over the log-resummed distri- 
bution with/without renormalon subtractions (red trian- 
gles and green squares). Both expansions are displayed 
including fixed order corrections up to order asijnz), 
a1{mz) and a'l{mz), as indicated by the orders 1, 2, 3, 
respectively. We immediately notice that the resummed 
result is not as effective in reproducing the total cross sec- 
tion as the fixed order expansion. Predictions that sum 
large logarithms have a substantial (perturbative) nor- 
malization uncertainty. On the other hand, as shown in 
Ref. j8], the resummation of logarithms combined with 
the profile function approach leads to a description of 
the thrust spectrum that converges nicely over the whole 
physical r range when the norm of the spectrum is di- 
vided out, a property not present in the spectrum of the 
fixed order expansion. 

In Fig. [2I the expansions of the partonic moments Mi , 
M2, and M3 are displayed in the fixed order expansion 
(blue circles) and the log-resummed result with cither 
the fixed order normalization (green squares) or a prop- 
erly normalized spectrum (red triangles). We observe 
that the fixed order expansion has rather small varia- 
tions from scale variation, but shows poor convergence 
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FIG. 3: Theory scan for uncertainties in pure QCD with massless quarks. The panels are fixed order (top- left), resummation 
without the nonperturbative correction (top-right), resummation with a nonperturbative function using the MS scheme for Q\ 
(bottom-left), resummation with renormalon subtraction and a nonperturbative function in the Rgap scheme for fli (bottom- 
right). 



indicating that its renormalization scale variation un- 
derestimates the perturbative uncertainty. For Mi the 
fixed order and log-resummed expressions with a com- 
mon fixed-order normahzation (blue circles and green 
squares) agree well at each order, indicating that, as ex- 
pected, large logarithms do not play a significant role for 
this moment. On the other hand, the expansion based on 
the properly normalized log-resummed spectrum exhibits 
excellent convergence, and also has larger perturbative 
uncertainties at the lowest order. In particular, for the 
red triangles the higher order results are always within 
the l-fj uncertainties of the previous order. The result 
shows that using the normalized log-resummed spectrum 
for thrust, which converges nicely for all r, also leads to 
better convergence properties of the moments. At third 
order all the fixed order and resummcd partonic moments 
are consistent with each other. Since the log-resummed 
moments exhibit more realistic estimates of perturbative 
uncertainties at each order, we will use the normalized 



resummed moments for our fit analysis.® 

In Fig. [3] we show how the inclusion of various ingredi- 
ents (fixed order contributions, log resummation, power 
corrections, renormalon subtraction) affects the conver- 
gence and uncertainty of our theoretical prediction for the 
first moment of the thrust distribution as a function of 
Q. From these plots we can observe four points: i) Fixed 
order perturbation theory does not converge very well, 
ii) Resummation of large logarithms in the distribution, 
when normalized with the integral of the resummed dis- 
tribution, improves convergence for every center of mass 
energy, iii) The inclusion of power corrections has the 



At N'^LL in our most complete theory set up the norm of the 
distribution and total hadronic cross section are fully compat- 
ible within uncertainties, so it does not matter which is used. 
Following Ref. fSj, at N^LL we choose to normalize the distribu- 
tion with the fixed-order total hadronic cross section since it is 
faster. 
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order as{mz) (with Q,i ) as{mz) (with fif^^^) 



NLL 

NNLL 
N^'LL (full) 



0.1173(82)(13) 
0.1159(41)(14) 
0.1153(21)(14) 



0.1172(82)(13) 
0.1139(15)(13) 
0.1140(07)(14) 



N LL(QCD + m6) 
N^LL(purc QCD) 



0.1160(20)(14) 
0.1156(21)(14) 



0.1146(07)(14) 
0.1142(07)(14) 



TABLE I: Central values for as{mz) at various orders with 
theory uncertainties from the parameter scan (first value in 
parentheses), and experimental and hadronic error added in 
quadrature (second value in parentheses). The bold N^LL 
value above the line is our final result, while values below 
the line show the effect of leaving out the QED and 6-mass 
corrections. 



FIG. 4: Difference between theoretical predictions with de- 
fault parameters for the first moment as function of Q when 
varying one parameter at a time. The red solid line corre- 
sponds to varying Aas(mz) = ±0.001 and the blue dashed 
lines to varying AS7i = ±0.1, with respect to the pure QCD 
best-fit values. There is a strong degeneracy of the two param- 
eters in the region Q > 100 GeV, which is obviously broken 
when considering values of Q below 70 GeV. 

effect of a 1/Q-modulated vertical shift on the value of 
the first moment, iv) The subtraction of the renormalon 
ambiguity reduces the theoretical uncertainty. This pic- 
ture for the first moment is consistent with the results of 
Ref. [8] for the thrust distribution. 

Another important element of our analysis is that we 
perform global fits, simultaneously using data at a wide 
range of center of mass energies Q. This is motivated by 
the fact that for each Q there is a complete degeneracy 
between changing as{mz) and changing fii, which can 
be lifted only through a global analysis. Fig. |4] shows 
the difference between the theoretical prediction of Mi 
as a function of Q, when as{mz) or ili are varied by 
±0.001 and ±0.1 GeV, respectively. We see that the 
effect of a variation in as{mz) can be compensated with 
an appropriate variation in f2i at a given center of mass 
energy (or in a small Q range) . This degeneracy is broken 
if we perform a global fit including the wide range of Q 
values shown in the figure. 

Finally, in Fig. [5] we show as{mz) extracted from fits 
to the first moment of the thrust distribution at three- 
loop accuracy including sequentially the different effects 
our code has implemented: 0(q!^) fixed order, N^LL re- 
summation, power corrections, renormalon subtraction, 
b-quark mass and QED. The error bars of the first two 
points at the left hand side do not contain an estimate 
of uncertainties associated with the power correction. 
Though smaller, the resummed result is compatible at 
the l-tJ level with the fixed order result. The inclusion of 
the power correction is the element which has the great- 
est impact on as{mz)', for the MS definition of ili it 
reduces the central value by 7%. The subtraction of the 



order Hi (MS) [GeV] fii (Rgap) [GeV] 

NLL 0.504(157)(45) 0.500(153)(45) 

NNLL 0.405(82) (47) 0.413(43) (44) 

N^LL (full) 0.318(75)(49) 0.377(39)(44) 

N3LL(qcd+™,) 0.310(74)(49) 0.369(34)(44) 

N^'LLcpurc QCD) 0.350(67) (49) 0.402(35) (44) 

TABLE IL Central values for Qi at the reference scales 
Ra = MA = 2 GeV and for Qi and at various orders. The 
parentheses show theory uncertainties from the parameter 
scan, and experimental and hadronic uncertainty added in 
quadrature, respectively. The bold value above the line is our 
final result, while the N^LL values below the horizontal line 
show the effect of leaving out the QED and fe-mass corrections. 







xV(dof) 


N^LL with fif'^^'P 


0.1140(07)(14) 


1.33 


N^LL with nf ^ 


0.1153(21)(14) 


1.33 


N'^LL no power corr. 


0.1236(39)(03) 


2.03 


0{a^) fixed order 


0.1305(39)(04) 


2.52 


no power corr. 







TABLE IIL Comparison of first moment fit results for anal- 
yses with full results and fii — Q^^^^, with Qi and no renor- 
malon subtractions, without power corrections, and at fixed 
order without power corrections or log resummation. The 
first number in parentheses corresponds to the theory uncer- 
tainty, whereas the second corresponds to the experimental 
and hadronic uncertainty added in quadrature for the first 
two rows, and experimental uncertainty for the last two rows. 



renormalon ambiguity in the Rgap scheme reduces the 
theoretical uncertainty by a factor of 3, while b-quark 
mass and QED effects give negligible contributions with 
current uncertainties. 
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cusimz) from global first moment thrust fits 



0.135- 
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0.125- 



0.120- 



0.115 



©(ttj) fixed-order 
0.1299 ±0.0038 



All errors: a,{mz) = 0.1140 ± 0.0016 



N^LL summation 
0.1245 ±0.0038 



perturbative error — 



+ Power Correction 
0.1156 ±0.0021 



perturbative error 




+ R- scheme 

0.1142 ± 0.0007 +b-mass&QED 
0.1140 ±0.0007 




FIG. 5: Evolution of the best-fit values for as{mz) from thrust first moment fits when including various levels of improvement 
with respect to fixed order QCD. Only points at the right of the vertical dashed line include nonperturbative effects. 
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FIG. 6: Distribution of best-fit points in the as(jnz)-1^\ and OLs(raz)-1^\ planes. The left panel shows results including 
perturbation theory, resummation of the logs, the soft nonperturbative function, and fii defined in the Rgap scheme with 
renormalon subtractions. The right panel shows the same results, but with VL\ defined in the MS scheme, and without 
renormalon subtractions. In both panels the dashed lines corresponds to an ellipse fit to the contour of the best-fit points 
to determine the theoretical uncertainty. The respective total (experimental -f theoretical) 39% CL standard error ellipses are 
displayed (solid lines), which correspond to l-cr (68% CL) for either one-dimensional projection. 



B. Uncertainty Analysis 



In Fig. [6] we stiow tlie result of our tlieory scan to 
determine the perturbative uncertainties. At each or- 



der we carried out 500 fits, with theory parameters ran- 
domly chosen in the ranges given in Table [VlII| of App. |X] 
(where further details may be found). The left panel of 
Fig. [6] shows results with renormalon subtractions using 
the Rgap scheme for 17 1, and the right-panel shows re- 
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FIG. 7: Experimental Ax'^ ~ 1 standard error ellipse (dotted 
green) at N^LL accuracy with renormalon subtractions, in 
the Qs-S S7i plane. The dashed blue ellipse represents the 
theory uncertainty which is obtained by fitting an ellipse to 
the contour of the distribution of the best-fit points. This 
ellipse should be interpreted as the I-ct theory uncertainty 
for 1-parameter (39% confidence for 2-parameters) . The solid 
red ellipse represents the total (combined experimental and 
perturbative) uncertainty ellipse. 



suits in the MS scheme without renormalon subtractions. 
Each point in the plot represents the result of a single 
fit. As described in App.|A] in order to estimate pertur- 
bative uircertainties, we fit an eUipse to the contour of 
best-fit points in the ols-1 f2i plane, and we interpret this 
as l-cr theoretical error ellipse. This is represented by 
the dashed lines in Fig. |6] The solid lines represent the 
combined (theoretical and experimental) standard error 
ellipses. These are obtained by adding the theoretical and 
experimental error matrices which determined the indi- 
vidual ellipses. The central values of the fits, collected in 
Tables |T] and [llj are determined from the average of the 
maximal and minimal values of the theory scan, and are 
very close to the central values obtained when running 
with our default parameters. The minimal values for 
these fits are quoted in Table |II1| as well. The best fit 
based on our full code has x^/dof = 1.325 ± 0.002 where 
the range incorporates the variation from the displayed 
scan points at N'^LL. The fit results show a substantial 
reduction of the theoretical uncertainties with increasing 
perturbative order. Removal of the ©(Aqcd) renormalon 
improves the perturbative convergence and leads to a re- 
duction of the theoretical uncertainties at the highest or- 
der by a factor of 2 in f2i, and factor of 3 in as{mz) 

To analyze in detail the experimental and the total un- 
certainties of our results, we refer now to Fig.[7j Here we 
show the error ellipses for our highest order fit, which in- 



cludes resummation, power corrections, renormalon sub- 
traction, QED and b-quark mass contributions. The 
green dotted, blue dashed, and the solid red lines rep- 
resent the standard error ellipses for, respectively, exper- 
imental, theoretical, and combined theoretical and ex- 
perimental uncertainties. The experimental and theory 
error ellipses are defined by Ax^ = 1 since we are most 
interested in the 1-dimensional projection onto as- The 
correlation matrix of the experimental, theory, and total 
error ellipses are (i,j = as, 2^1) 



T/thco 



1.93(15) -10" 

l-lQ-^'GcV 0.79(13) ■ 

1.85 • 10-5 GeV 



-1.18(13) • 10"'' GeV \ 
GeV 0.79(13) • lO^^ GeV^ ) ' 



-1.18(13) 

5.56-10-^ _ 
1.85 • 10-5 GeV 5.82 • lO^^ QeV^ 

1-10-6 -0.99(13) • lO-'^ GeV ' 



2.49(15) • lu — u.yai^io; • lu vjc v 

-0.99(13) • 10--* GeV 1.37(13) • IQ-^ GcV^ 



where the experimental correlation coefficient is signifi- 
cant and reads 



Pan 



0.96(14) . 



(30) 



Adding the theory scan uncertainties reduces the corre- 
lation coefficient in Eq. ( 30 ) to 



= -0.54(8). 



(31) 



In both Eqs. (30) and (31 1 the numbers in parentheses 



capture the range of valu es o btained from the theory 



scan. From V^^^ in Eq. (29) it is possible to extract 



the experimental uncertainty ior as and f2i and the un- 
certainty due to variations of f2i and a^, respectively: 



(Ta.xH-pin^ 0.0004. 



(32) 



C'^= ^0,^1- ^'0= 0.013 GeV, 

^al = I Pan I = 0.0014, 
't'^, = Wn\ = 0.044 GeV. 

Fig. [7] shows the total uncertainty in our final result 
quoted in Eq. (34) below. 



The correlation exhibited by the green dotted experi- 
mental error ellipse in Fig.[7]is given by the line describing 
the semimajor axis 



32.82 GeV 



0.1255 -as(mz). 



(33) 



Note that extrapolating this correlation to the extreme 
case where we neglect the nonperturbative corrections 
(r^i = 0) gives as{mz) 0.1255. 
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C. Effects of QED and the &-mass 

The experimental correction procedures applied to the 
AMY, JADE, SLC, DELPHI and OPAL data sets were 
typically designed to eliminate initial state photon radi- 
ation, while those of the TASSO, L3 and ALEPH coUab- 
orations eliminated initial and final state photon radia- 
tion. It is straightforward to test for the effect of these 
differences in the fits by using our theory code with QED 
effects turned on or off depending on the data set. Using 
our N'^LL order code in the Rgap scheme we obtain the 
central values as{mz) — 0.1143 and ili = 0.376 GeV. 
Comparing to our default results given in Tabs. P and [Til 
which are based on the theory code were QED effects are 
included for all data sets, we see that the central value 
for as is larger by 0.0003 and the one for f2i is smaller 
by 0.001 GeV. This shift is substantially smaller than 
our perturbative uncertainty. Hence our choice to use 
the theory code with QED effects included everywhere 
as the default for our analysis does not cause an observ- 
able bias regarding experiments which remove final state 
photons. 

By comparing the N^LL (pure massless QCD) and 
N'^LL (QCD + TTif,) entries in Tabs.|l]and[ll]we see that in- 
cluding finite &-mass corrections causes a very mild shift 
of ~ +0.0004 to as{mz), and a somewhat larger shift 
of ~ —0.033 GeV to Jli. In both cases these shifts are 
within the 1-a theory uncertainties. In the N'^LL (pure 
massless QCD) analysis the &-quark is treated as a mass- 
less flavor, hence this analysis differs from that done by 
JADE |23| where primary b quarks were removed using 
MC generators. 



D. Final Results 

As our final result for as{mz) and fJi, obtained at 
N'^LL order in the Rgap scheme for ^li(/iA, A^a), includ- 
ing bottom quark mass and QED corrections we obtain 

as{mz) = 0.1140 ± (0.0004)exp (34) 

± (0.0013)hadr ± (0.0007)pert, 

f^i(i?A,MA) = 0.377 ± (0.013)cxp 

± (0.042)„^(^,) ± (0.039)pert GeV, 

where i?A = Ma = 2 GeV and we quote individual 1-a 
uncertainties for each parameter. Here x^/dof = 1.33. 
Eq. pi] ) is the main result of this work. 

In Fig. [8] we show the first moment of the thrust dis- 
tribution as a function of the center of mass energy Q, 
including QED and m?, corrections. We use here the best- 
fit values given in Eq. (34). The band displays the theo- 



Mi(Q) 



T 1 1 1 1 1 1 1 1- 

Fit at N^^LL for Q'j((7!z) and Hi 
theory scan error 




150 200 
Q (GeV) 



FIG. 8: First moment of the thrust distribution as a func- 
tion of the center of mass energy Q, using the best-fit values 
for asijnz) and fii in the Rgap scheme as given in Eq. (34 1. 
The blue band represents the perturbative uncertainty deter- 
mined by our theory scan. Data is from ALEPH, OPAL, L3, 
DELPHI, JADE, AMY and TASSO. 




retical uncertainty and has been determined with a scan 
on the parameters included in our theory, as explained in 
App. [X] The fit result is shown in comparison with data 



0.II2 0.II3 0.II4 0.II5 0.II6 0.II7 

FIG. 9: Comparison of as{mz) and Sli determinations from 
thrust first moment data (red upper right ellipses) and thrust 
tail data (blue lower left ellipses). The plot corresponds to 
fits with N'^LL accuracy and in the Rgap scheme. The tail 
fits are performed with our improved code which uses a new 
nonsingular two-loop function, and the now known two-loop 
soft function. Dashed lines correspond to theory uncertain- 
ties, solid lines correspond to A.x^ — 1 combined theoretical 
and experimental error ellipses, and wide-dashed lines corre- 
spond to Ax'^ ~ 2.3 combined error ellipses (corresponding 
to 1-a uncertainty in two dimensions). 



from ALEPH, OPAL, L3, DELPHI, JADE, AMY and 
TASSO. Good agreement is observed for all Q values. 
It is interesting to compare the result of this analysis 



13 



with the result of our earher fit of thrust tail distributions 
in Ref . [8] . This is shown in Fig. |9] Here the red upper 
shaded area and corresponding ellipses show the results 
from fits to the first moment of the thrust distribution, 
while the blue lower shaded area and ellipses show the 
result from fits of its tail region. Both analyses show 
the theory (dashed lines) and combined theoretical and 
experimental (solid lines) standard error ellipses, as well 
as the ellipses which correspond to Ax^ = 2.3 (68% CL 
for a two-parameter fit, wide-dashed lines). We see that 
the two analyses are compatible. 



order 0{a'i) 0(a?) 

(i) Rgap R-RGE 0.1159(27)(14) 0.1146(06)(14) 

(ii) Rgap FO Subt. 0.1185(63)(15) 0.1138(20)(14) 
(ill) MS for Hi 0.1278(124)(19) 0.1186(38)(14) 

TABLE IV: MS scheme values for as{mz) obtained from var- 
ious fixed order analyses. The first value in parentheses is the 
uncertainty from higher order perturbative corrections (ob- 
tained by the method described in the text), while the second 
value is the combined experimental and hadronization uncer- 
tainty. 



IV. FIXED ORDER ANALYSIS OF Mi 

It is interesting to compare the result of our best fit 
with an analysis where we do not perform resummation 
in the thrust distribution, but where power corrections 
and renormalon subtractions are still considered. This 
is achieved by setting the scales ^ih, MSj MJi /^ns hi^ our 
theoretical prediction all to a common scale /i ^ Q. We 
use R for the scale of the renormalon subtractions and 
renormalization group evolved power correction. Finally 
we will neglect QED and 6-mass corrections in this sub- 
section. Up to the treatment of power corrections and 
perturbative subtractions, the fixed order results used for 
this analysis are thus equivalent to those used in Ref. [JS] . 

The OPE formula for the first moment in the Rgap 
scheme for this situation is given by 

Ah = Mfs^P (i?, ^l) + ^IM^ , (35) 
ni{R, A{R, fi) - A(i?A, ma) , 



order 



o(«? 



In Eq. (35 1, the Qi with no arguments is the value de- 



termined by the fits, which is in the Rgap scheme at the 
reference scale = Ra = 2 GeV. Here A{R, n) is the 
running gap parameter, and A(i?, /z) — A(i?A, A*a) is used 
to sum logarithms from (i?A,MA) to (i?, m) in Eq. (35 1. 
The analytic expression for A{R, fj,) — A(f?A, Ma) can be 
found in Eq. (41) of Ref. [S] (see also HH). The per- 
turbative M^^'^P is related to the perturbative MS result 
by 



Mfs-P(i?,Ai) = Mf«(/^) + 



MS/ 



2d{R,ti) 



(36) 



where the subtractions terms are [H [12] 

Si (R, m) = -0.848826Lfl , (37) 
S2iR,fi) = -0.156279 - 0.46663Lfl- 0.517864L^, 
SsiR, fJ-)^- 0.552986 - 0.622467Li^ - 0.777219L|j 



- 0.421261L' 



R ' 



with Lji = ln(/i/i?). In Eq. (|36j) 6{R, ^) cancels the 
©(Aqcd) renormalon in M^^(/i), and it is crucial that 



(i) Rgap R-RGE 0.407(8)(45) 0.400(8)(45) 
(fi) Rgap FO Subt. 0.216(126)(133) 0.359(42)(62) 
(iii) MS for Hi 0.388(62) (47) 0.350(54) (44) 

TABLE V: Qi or fli values obtained from fixed order analyses 
at various orders. The first value in parentheses is the uncer- 
tainty from higher order perturbative corrections (obtained 
by the method described in the text), while the second value 
is the combined experimental and hadronization uncertainty. 



the coupling expansions in both these objects are done 
at the same scale, as{fJ,), for this cancellation to take 
place. The relation to the MS scheme power correction 
is r^i = rii -I- 5(i?A, Ma), and the OPE in the MS scheme 
at this level is 



Ml =Mf^(Ai) + 



2rji 



(38) 



In the MS result there are no perturbative renormalon 
subtractions (and thus no log resummation related to 
the renormalon subtractions) and the parameter Qi has 
a Aqcd renormalon ambiguity. 

We will perform fits to the experimental data follow- 
ing the same procedure discussed in the previous section. 
Using Eq. (35 1 we consider two cases, i) R ^ Q where fii 
is renormalization group evolved to R and there are no 
large logarithms in the renormalon subtractions, and ii) 
fixing R at the reference scale, R — 2 GeV, in which case 
large logarithms are present in the renormalon subtrac- 
tions. We will also consider a third case, iii), using the 
MS-OPE of Eq. J38| . Results for these fits are shown 
in Tabs. [Tvl and [VpFor ah cases xVdof ~ 1.32. 

Q, so there are no 



For case i) we take R 
large logarithms in the 5{R,fi) of Eq. (35), and all large 



M 



logarithms associated with renormalon subtractions are 
summed in A(i?, /i) — A(i?A,MA)- Here we estimate 
the perturbative uncertainty in as{mz) and fli by vary- 
ing the renormalization scale fj, and the scale R inde- 
pendently in the range {2Q,Q/2}. We use one-half 
the maximum minus minimum variation as the uncer- 
tainty, and the average for the central value. The re- 
sults for both as{mz) and J7i are fully compatible at 
1-(T to our final results shown in Eq. (34). The agree- 
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ment is even closer to the central values for the fits with- 
out QED or 6-mass corrections in Tabs. [T| and [iTl namely 
cxsimz) = 0.1142(07)(14) and Oi = 0.402(35)(44). The 
one difference is that the perturbative uncertainty for Qi 
m Tab.|v] is a factor of three smaller. The case i) results 
in the table also exhibit nice order-by-order convergence, 
and if one plots Mi versus Q (analogous to Fig. [2]) the 
uncertainty bands are entirely contained within one an- 
other. In order to be conservative, we take our resum- 



mation analysis in Eq. (34) as our final results (with its 



larger perturbative uncertainty and inclusion of QED and 
&-mass corrections). 

For case ii) we take i? ~ 2 GeV and fi ^ Q a,s typi- 
cal values, so there are large logarithms, ln{R/Q), in the 
S{R, fi) renormalon subtractions. The central value for 

(mz) at ©(a^) is again fully compatible with that in 
Eq. ( |34[ ). Here we estimate the perturbative uncertainty 
in asljnz) by varying /i G {2 Q, Q/2} and R ~ 2±1 GeV. 
Due to the large logarithms the perturbative uncertainty 
in ag(mz) for case ii), shown in Tab. IV is three times 
larger than for case i) . It is also compatible with the dif- 
ference between central values at 0{al) and 0{al). To 
estimate the uncertainty for fli we only vary which 
leads to the rather large error estimate for fii shown in 
Tab.fVl The contrast between the precision of the results 
in case i), to the results in case ii), illustrates the im- 
portance of summing large logarithms in the renormalon 
subtractions. 

For case iii), where the f2i power correction is defined 
in MS we do not have renormalon subtractions (and 
hence no large logs in subtractions). Due to the poor 
convergence of the fixed order prediction for the first mo- 
ment, seen from the blue fixed order points in Fig. [2] it is 
not clear whether varying fi in the range {2 Q, Q/2} gives 
a realistic perturbative uncertainty estimate. Hence we 
determine the perturbative uncertainty for case iii) in 
Tabs. IV and [v] by varying ^ in the range {2(5, Q/2} 
and multiply the result by a factor of two. The pertur- 
bative uncertainties for asijnz) are a factor of two larger 
than in case ii). The central values for as{mz) in case 
iii) are also larger, but are compatible with those in case 
ii) and Eq. (34) within l-cr. 

It is interesting to compare our results to those of 
Ref. j48], which also performs a fixed order analysis 
at 0{al), and incorporates subtractions based on the 
dispersive model. ^ Here the subtractions contain loga- 
rithms, \n(fii/ fj,), where fij ^ 2 GeV and fi ^ that are 
not resummed. From a fit to Mi in thrust they obtained 
as{mz) = 0.1166 ± 0.0015exp ± 0.0032th where the first 
uncertainty is experimental and the second is theoreti- 



® On the experimental side, Ref. |48l uses only the new JADE data 
from I23| and OPAL data. In our analysis the new JADE was 
excluded, but we utilized a larger dataset that includes ALEPH, 
OPAL, L3, DELPHI, AMY, TASSO, and older JADE data. This 
may have a non-negligible impact on the outcome of the com- 
parison. 



cal. Our corresponding result is the one in case ii), and 
the central values and uncertainties for as{mz) are fully 
compatible. The perturbative uncertainty they obtain is 
a factor of 1.6 larger than ours. It arises from varying the 
renormalization scale ^ E {2Q,Q/2}, the 0{a^) Milan 
factor A4 by 20%, and the infrared scale /i/ = 2 ± 1 GeV 
in the dispersive model. In our analysis there is no pre- 
cise analog of the Milan factor because our subtractions 
and Rgap scheme for fii fully account for two and three 
gluon infrared effects up to 0{al) that are associated to 
thrust. Other than this, the difference can be simply at- 
tributed to the differences in subtraction schemes which 
have an impact on the /i scale uncertainty. Finally, note 
that we have implemented the analytic results of Ref. [48] 
and confirmed their fi and fij uncertainties. 



V. JADE DATASETS 

As discussed in Sec. [l] our global dataset includes 
thrust moment results from ALEPH, OPAL, L3, DEL- 
PHI, AMY, TASSO and the JADE data from Ref. ^22]. 
In this section we discuss the impact on the results 
in Sees. 



Ill and IV of replacing the JADE data from 



Ref. 



with moment results from an updated analysis 
carried out in Ref. |23| . which removes the contributions 
from primary bb pair production and provides in addi- 
tion measurements at Q = 14 and 22 GeV. In Fig 



10 



we 

show the data for Mi, including the JADE results from 
Refs. [H] and . The most significant difference occurs 
at Q = 44 GeV. Our analysis will treat these datasets on 
the same footing without attempting to account for the 
effect of removing the 66' s. 

For our analysis here, with theory results at 
N'^LL -|- 0{al), we continue to exclude center of mass en- 
ergies Q < 22 GeV as in Sec. |III| The dependence of the 
global fit result on the data set for Mi is shown in Fig. 11 
Theoretical uncertainties are analyzed again by the scan 
method giving the central dots and three inner ellipses, 
while the outer three ellipses show the respective com- 
bined l-cr total experimental and theoretical uncertain- 
ties. Using all experimental data but excluding JADE 
measurements entirely gives the fit result shown by the 
upper blue ellipse. This result is compatible at l-cr with 
the central red ellipse which shows our default analysis, 
using the Ref. [22^ JADE Mi measurements. Replacing 
these two JADE data points by the four Q > 22 GeV 
JADE Ml results from Ref. yields the lower green 
ellipse (whose center is ~ 1.5-cr from the central ellipse). 
For this fit the x^/dof increases from 1.33 to 1.52 demon- 
strating that there is less compatibility between the data. 
For this reason, together with the concern about the im- 
pact of removing primary 66 events with MC simulations, 
we have used only JADE data from Ref. [22] in our main 
analysis. 

A similar pattern is observed using the fixed order fits 
of Ml discussed in Sec. |IV| In this case it is also straight- 
forward to include the Q = 14, 22 GeV JADE data from 
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FIG. 10: Experimental data for the first moment of thrust. 
The solid line corresponds to the result from the first row of 
Tab. |IV[ and uses a fixed order code with power corrections in 
a renormalon-free scheme, but no resummation (neither QED 
nor bottom mass corrections). 



Ref. 23_. If these two points are added to our default 
dataset (which contains Q = 35 and 45 GeV as the lowest 
g results for Mi) then we find a^imz) = 0.1155 ±0.0012 
and Qi = 0.361 ± 0.035 GeV with xVdof = 1.3. This 
is compatible at 1-cr with our final pure QCD result in 
Tab. m If we include the entire set of JADE data from 
Ref. [23] instead of those from Ref. [22] then we find 
as{mz) = 0.1166 ± 0.0012 and = 0.306 ± 0.033 GeV 
with x^/dof = 1.6, very similar to the values observed 

Hence, overall the 



for the green lower ellipse in Fig. 11 



fixed order analysis does not change the comparison of 
fits with the two different JADE datasets. 



VI. HIGHER MOMENT ANALYSIS 

In this section we consider higher moments, Af„>2, 
which have been measured experimentally up to n = 5. 
From Eq. (21) we see that these moments have power 
corrections oc 1/Q'^ for A: > 1. Since for the perturbative 
moments we have Mn/Mn+i — 4-9, we estimate that 
the 1 /Q^ power corrections are suppressed by 9Aqcd /Q 
which varies from 1/8 to 1/44 for the Q- values in our 
dataset, Q > 35 GeV. Hence, for the analysis in this 
section we can safely drop the l/Q^ and higher power 
corrections and use the form 



2nni 



(39) 



By using our fit results for as{mz) and fii from 



Eq. (34) we can directly make predictions for the mo- 
ments M2_3_4^5. This tests how well the theory does at 
calculating the perturbative contributions M2,3,4,5. The 
results for these moments are shown in Fig. [12] and cor- 
respond to xVdof = 1.3, 2.5, 0.8, 1.1 for n = 2, 3, 4, 5 re- 
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FIG. 11: Fit results when using ALEPH, DELPHI, OPAL, L3, 
AMY, TASSO, but no JADE data (upper blue ellipse), when 
also including JADE data from Ref. [22) (red central ellipse) 
[ our default data set ] , and when instead including the JADE 
data from Ref. (green lower ellipse). The ellipses here 
correspond to 1-a for two parameters (68% CL). 



spectively, indicating that our formalism does quite well 
at reproducing these moments. The larger x^/dof for 
n = 3 is related to a quite significant spread in the ex- 
perimental data for this moment at Q > 190 GeV. Note 
that we also see that the relation Af„/A/„+i ~ 4-9 is 
satisfied by the experimental moments. 

An alternate way to test the higher moments is to 
perform a fit to this data. Since we have excluded the 
new JADE data in Ref. [53], we do not have a signifi- 
cant dataset at smaller Q values for the higher moments. 
With our higher moment dataset the degeneracy between 
as{mz) and fii is not broken for n > 2, and one finds 
very large experimental errors for a two-parameter fit al- 
ready at n = 2. However we can still fit for as{mz) from 
data for each individual Mn>2 by fixing the value of fli 



to the best fit value in Eq. (34) from our fit to Mi. For 
this exercise we use our full N^L + 0{a^) code, but with 
QED and mass effects turned off. The outcome is shown 



in Fig. 13 and Tab. VI We find only a little dependence 
of ag on n, and all values are compatible with the fit to 
the first moment within less than 1-cr. This again con- 
firms that our value for Qi and perturbative predictions 
for Mn>2 are consistent with the higher moment data. 

In Ref. (48] a two-parameter fit to higher thrust mo- 
ments was carried out using OPAL data and the latest 
low energy JADE data. For n = 2 to n = 5 the results 
increase linearly from as{mz) = 0.1202 ± (0.0018)oxp ± 
(0.0046)th to as{mz) 0.1294 ± (0.0027)oxp ± (0.0070)th 
respectively, and the weighted average for the first five 



moments of thrust is as{mz) 



0.1208 ± 0.0018cxp ± 
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FIG. 12: Predictions for the higher moments M2, M3, M4, 
M5 using the best fit values from Eq. (341, and our full 
N^LL + 0(0:3) code in the Rgap scheme, but with QED and 
mass effects turned off. The central points use different sym- 
bols for different moments. 



n 


as{mz) 


Ath[as] 


Acxp[as] 


xVdof 


2 


0.1149 


0.0009 


0.0005 


1.24 


3 


0.1157 


0.0009 


0.0005 


1.87 


4 


0.1151 


0.0011 


0.0010 


0.39 


5 


0.1156 


0.0015 


0.0010 


0.23 



TABLE VI: Numerical results for Us from one-parameter fits 
to the Mn moments. The second column gives the central 
values for as{mz), the third and fourth show the theoretical 
and experimental errors, respectively. Since Q,i was fixed for 
this analysis we do not quote a hadronization error. 



0.0045th. The results are fully compatible within the 
uncertainties, and there is an indication of a trend to- 
wards larger as{mz) extracted from higher moments. In 
our analysis we do not observe this trend, but our re- 
sults should not be directly compared since we have only 
performed a one parameter fit. After further averag- 
ing over results obtained from event shapes other than 
thrust Ref. [IS] obtained as their final result as{mz) = 
0.1153 ± 0.0017oxp ± 0.0023th. This is again perfectly 



compatible with our result in Eq. (34) 



5- 



4- 



3- 



2- 



1 - 



0.110 0.112 0.114 0.116 0.118 0.120 

FIG. 13: One-parameter fits for as{mz) to the first five mo- 
ments. We use our full set up with power corrections and 
renormalon subtractions, but with QED and mass corrections 
turned off. The value of Oi is fixed from Eq. ( 34 1 . The er- 



ror bars include theoretical and experimental errors added in 
quadrature (not including uncertainty in ili). 



rections and their constraints from experimental data. 
There are two types of power corrections that are relevant 
for the cumulants, those defined rigorously by QCD ma- 
trix elements which come from the leading thrust factor- 
ization theorem, f2^, and those from our simple parame- 
terization of higher order power corrections in Eq. (|15[) , 



^n,j>i- For the latter a systematic matching onto QCD 
matrix elements has not been carried out and the corre- 
sponding perturbative coefficients have not been deter- 
mined. 

For the second cumulant AI2 both types of power cor- 
rection contribute to the leading term in the com- 
bination 



Q'2 = n'2 + Ml 1 fii 1 , 



(40) 



Without a calculation of the perturbative coefficient Af 1.1 
we cannot argue that either one dominates, and hence we 
keep both of them. In terms of this parameter the OPE 
with its leading power correction for the second cumulant 
becomes simply 



VII. HIGHER POWER CORRECTIONS FROM 
CUMULANT MOMENTS 



(41) 



In this section we use cumulant moments as defined 



in Eq. (27) to discuss the presence of higher power cor- 



where M2 computed from our leading order factoriza- 
tion theorem, see Eq. (11) 



For the third cumulant Afg 
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FIG. 14: Prediction of cumulants using our best-fit values 
for as{mz) and Qi from the fit to the first thrust moment. 
The band includes only the theoretical uncertainty from the 
random scan. The theory prediction includes QED and mass 
corrections, in contrast to our numerical analysis which has 
no QED and 6-mass effects and uses our default model, which 
translates into the following values for higher nonperturbative 
power corrections: Q'2 = ^1/4, fig = Qf/S, fii = 3Qt/32, 
n's = 3f2f/32. 



the power correction from the leading thrust factoriza- 
tion theorem is while that from the subleading 



factorization theorem is l/Q , so 



6M2,il7i 



8^ 



(42) 



where we keep both of these power corrections. 

For our analysis we assume that the perturbative coef- 
ficients Afi.i and M2,i get contributions at tree-level, 
and hence that their logarithmic dependence on Q is 
Q!s-suppressed. Thus for fits to M2 and A/3 we con- 
sider the three parameters ^'2, M2,i^i,i, and £^3. Our 
theoretical expectations are that ($7^)^/" ~ Aqcd and 

Since most of the experimental collaborations pro- 
vide measurements only for moments we computed the 
cumulants using Eq. ^. To propagate the errors to 
the ri-th cumulant one needs the correlations between 
the first n moments, both statistical and systemati- 
cal. Following experimental procedures we estimate the 
statistical correlation matrix from Monte Carlo simu- 
lations. These matrices are provided in Ref. |65j for 





central 


Ath 


^cxp 


X 

dof 




0.74 


0.09 


0.11 


0.72 


(62) ^/^ 

(63) ^/^ 


1.21 
-2.61 


0.10 
0.15 


0.22 
1.51 


0.93 



TABLE VII: Determination of power corrections from fits to 
M2 and A/3. All values in the table are in GeV. Columns two 
to four correspond to central value, theoretical uncertainty, 
and experimental uncertainty, respectively (the latter includes 
both statistical and systematic errors added in quadrature). 
The values displayed correspond to the linear combinations 
in Eq. (43 1, which for A/3 diagonalize the experimental error 
matrix. 



Q = 14, 91.3, 206.6 GeV.i° The computation of these 
matrices does not depend on the simulation of the de- 
tector and hence can be a priory employed on the data 
provided by any experimental collaboration. It was found 
that statistical correlation matrices depend very mildly 
on the center of mass energy, and our approach is to use 
the matrix computed at 14 GeV for Q < 60 GeV, the 
one computed at 91.3 for 60 GeV < Q < 120 GeV and 
the one at 206.6 GeV for Q > 120 GeV. The systematic 
correlation matrix for the moments is estimated using 
the minimal overlap model based on the systematic un- 
certainties, and then converted to uncertainties for the 
cumulants. We use this method even for the few cases in 
which experimental collaborations provide uncertainties 
for the cumulants directly, since we want to treat all data 
on the same footing. In these cases we have checked that 
the results are very similar. 

To some extent the prescription we employ lies in 
between two extreme situations: a) moments are com- 
pletely uncorrelated, and b) cumulants are completely 
uncorrelated. Situation a) corresponds to the naive as- 
sumption that the moments are independent. Situation 
b) is motivated by considering that properties like the lo- 
cation of the peak of the distribution (^ Afi), the width 
of the peak (~ A/j), etc. are independent pieces of in- 
formation. By assuming moments are uncorrelated one 
overestimates the errors of the cumulants. This would 
translate into larger experimental errors for our fit re- 
sults and very small x^/dof. Assuming that cumulants 
are uncorrelated induces very strong positive correlations 
between moments, which then leads to small uncertain- 
ties for the cumulants, especially for the variance, and 
larger x^/dof values. With the adopted prescription we 
use one finds a weaker positive correlation among mo- 
ments, which translates into a situation between these 
two extremes. 



We thank Christoph Pahl for providing details on the use of 
correlation matrices for moments. 

One might also construct the correlation matrices using the sta- 
tistical and systematic errors from the thrust distributions them- 
selves. Bins in distributions are statistically independent and 
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FIG. 15: Determination of power corrections from fits to data. On the vertical axes we display the n-th experimental cumulant 
with the perturbative part subtracted — M'„ . The error bars shown are experimental (statistical and systematic combined) 
added in quadrature with perturbative errors from the random scan over the profile parameters. The top-left panel shows the 
fit to 0,2/0^ , and the top-right panel shows the fit to M2,i and Q.'^/Q'^ through the linear combinations in 62,3. The 

bottom two panels for n = 4, 5 show a simple fit to Ma^Sli^i and A/4jf2i,i taking ^4 = ^5 = 0. 



For our analysis we use our highest order code as de- 
scribed in Sec. Ill and take the value as{mz) — 0.1142 



obtained in our fit to the first moment data with this code 
(see Tab.lTj). Since we are analyzing cumulants M4>2 the 
value of Hi is not required, and there is no distinction be- 
tween having this parameter in MS or the Rgap scheme. 
Hence in order to fit for higher power corrections we use 
our purely perturbative code in the MS scheme. Thus 
all of the power correction parameters extracted in this 
section are in the MS scheme. The perturbative error is 
estimated as in Sec. |III[ by a 500 point scan of theory 
parameters (see App. |A[). 

Before we fit for the higher power corrections, we will 
check how well our factorization theorem predicts the ex- 
perimental cumulants using a simple exponential model 



systematic correlations are estimated using the minimal overlap 
model. Unfortunately this can introduce biases, and we thank 
Christoph Pahl for clarifying this point. 



for the nonperturbative soft function (the model with 
only one coefficient cq = 1 from Refs. [51 |SSj). This 
model has higher power corrections that are determined 
by its one parameter fli: fl'2 = = 

= 317f/32, fl'r, = 3nl/32. Results are shown in 
Fig. [T4j where good agreement between theory and data 
is observed. 

For the in Fig. 
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we also observe that M'^.i/M^ 



1/10, so the (n -|- l)-tn order cumulant is generically one 
order of magnitude smaller than the n-th order cumulant. 

Next we will fit for the power correction parameters 
^'21 M2.1 ^1.1, and rig. For this analysis we neglect QED 
and 6-mass effects. To facilitate this we consider the dif- 
ference between the experimental cumulants and the 
perturbative theoretical cumulants M^, namely — M2 



and M3 — ML From Eqs. (41 1 and (42) these differences 



are determined entirely by the power correction parame- 
ters we wish to fit. The results are shown in Tab. lVIll and 
the upper two panels of Fig. 15 From the M2 — M2 fit 
a fairly precise result is obtained for (jlj)^^^- Its central 
value of 740 MeV is compatible with ~ 2Aqcd , and hence 
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agrees with naive dimensional analysis. Interestingly, we 
have checked that including a constant and l/Q term in 
the second cumulant fit one finds that their coefficients 
are compatible with zero, in support of the theoretically 
expected l/Q^-depcndcncc. 

For the fit to il/g — A/3 there is a strong correlation 
between fJg and Af2,i f^i.i even though they occur at 
different orders in 1/Q. Since the is quadratic in these 
two parameters we can determine the linear combinations 
that exactly diagonalize their correlation matrix: 



6, 



6Af2,i 



0.07 



(0.3105 GeV"^)r23 



3 ' 



(43) 



63 = - (0.3105 GeV) 



6Af2,i 
0.07 



1,1 



Note that these combinations arise solely from experi- 
mental data. We have presented the coefficients of these 
combinations grouping together a factor of (6Af2,i/0.07), 
which is close to unity if 6A/2 1 — Mi . The results in 



Tab. VII exhibit a reasonable uncertainty for O2, but a 

Hence, at this time it is not 

and 

A'/2,i independently. As in the previous case, the 
fit docs not exhibit any evidence for a 1/Q correction, 
confirming the theoretical prediction for this cumulant. 



large uncertainty for 63 
possible to determine the original parameters fig 



In Fig. 15 we also show results for cumulant differ- 
ences M'^ — versus Q ioi n — 4 and n = 5. In all 
cases n = 2, 3, 4, 5 the perturbative cumulants are the 
largest component of the cumulant moments M^, as can 
be verified by the reduction of the values by a factor of 
2-3 in Fig. [15] compared to the values in Fig.[T4] We also 
observe an order of magnitude suppression between the 
{n + iyth andn'th terms, {M^_^_l- M^_^_l)/{M^- M^) - 
1/10. For n = 4, 5 the OPE formula in Eq. ([27| involves 
both 2^^l'^/Q" terms and terms with non-trivial pertur- 
bative coefficients: {2nMn-i,i^i,i)/Q'^ + - ■ ■ (where here 
the ellipses are terms at and beyond). If the former 

dominated we would expect a suppression by 2 Aqcd /Q 
for the {n + l)'th versus n'th term. The observed sup- 
pression by 1/10 is less strong and is instead consistent 
with domination by the power correction terms in 

the n — A,b cumulant differences. This would imply 
[{n+l)Mn^i\/[nMn-i,i\ ~ 1/10 and could in principle be 
verified by an explicit computation of these coefficients. 
In Fig. [15^ we show fits to a 1/Q^ power correction, which 
are essentially dominated by the lowest energy point at 

the Z-pole. The results are JsM^^i 0.20 ± 0.08 



from fits to Af^ and ^10 Af4,i = 0.07±0.06 from fits 
to A/5. These values agree with our expectation of the 
~ 1/10 suppression between the two Af„.i perturbative 
coefficients. 

In this section we have determined the power 
correction parameter fig with 25% accuracy, and find 
it is 3.8 a different from zero. For the higher moments 
there are important contributions from a fii.i/Q^ power 



correction, which appears to even dominate for n > 4. 
Clearly experimental data supports the pattern expected 
from the OPE relation in Eq. E7\. 



VIII. CONCLUSIONS 

In this work we have used a full r-distribution factor- 
ization formula developed by the authors in a previous 
publication |8| to study moments and cumulant moments 
(cumulants) of the thrust distribution. Perturbatively 
it incorporates 0{a^) matrix elements and nonsingular 
terms, a resummation of large logarithms. In ' t, to N'^LL 
accuracy, and the leading QED and bottom mass correc- 
tions. It also describes the dominant nonperturbative 
corrections, is free of the leading renormalon ambiguity, 
and sums up large logs appearing in perturbative renor- 
malon subtractions. 

Theoretically there are no large logs in the perturbative 
expression of the thrust moments, and when normalized 
in the same way the perturbative result from the full t 
code with resummation agrees very well with the fixed or- 
der results. Nevertheless, when the code is properly self 
normalized it significantly improves the order-by-order 
perturbative convergence towards the 0{a^) result. In 
particular, the results remain within the perturbative er- 
ror band of the previous order, in contrast to what is 
observed using fixed order expressions. This lends sup- 
port to the theoretical uncertainty analysis from the code 
with resummation. 

From fits to the first moment of the thrust distribu- 
tion. Ml, we find the results for as{mz) and the leading 
power correction parameter f2i given in Eq. (34). They 



are in nice agreement with values from the fit to the tail 
of the thrust distribution in Ref. . The moment results 
have larger experimental uncertainties, and these domi- 
nate over theoretical uncertainties, in contrast with the 
situation in the tail region analysis of Ref. 8 . Repeating 
the Ml fit using a fixed order code with no Inr resum- 
mation, but still retaining the summation of large logs 
in the perturbative renormalon subtractions, yields fully 
compatible results for as{mz) and fli. 

Using a Fourier space operator product expansion we 
have parameterized higher order power corrections which 
are beyond the leading factorization formula, and an- 
alyzed the OPE both for moments Af„ and cumulants 
Af^. In the moments A/„ the Vti/Q power correction from 
the leading factorization theorem enters with a perturba- 
tive suppression in its coefficient, and dominates numeri- 
cally over higher 1/Q corrections. In contrast, the cumu- 
lants Af^>2 depend on higher order cumulant power cor- 
rections Cl'^/Q" from the leading factorization theorem, 
and are independent of ili/Q, fiJ^.^/Q"^^. Data 
on these cumulants appear to indicate that they receive 
important contributions from a 1/Q^ power correction 
that enters at a level beyond the leading thrust factor- 
ization theorem. Thus the OPE reveals that cumulants 
are appealing quantities for exploring subleading power 
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parameter default value range of values 





2GeV 


1.5 to 2.5 GeV 




5 


2 to 8 


t2 


0.25 


0.20 to 0.30 


ej 





-1, 0, 1 




1 


0.5 to 2.0 


Us 





-1, 0, 1 


T-icusp 
^ 3 


1553.06 


-1553.06 to +4659.18 


33 





-3000 to +3000 


S3 





-500 to +500 


£2 





-1, 0, 1 


£3 





-1, 0, 1 



TABLE VIII: Theory parameters relevant for estimating the 
theory uncertainty, their default values and range of values 
used for the theory scan during the fit procedure. 



corrections. We performed a fit to the second cumulant 
and determined a non- vanishing Cl^/Q^ power correction 
with a precision of 25%. 

It would be interesting to extend the analysis per- 
formed here, based on OPE formulas related to factoriza- 
tion theorems, to other event shape moments and cumu- 
lants. Examples of interest include the heavy jet mass 
event shape [3 angularities [701 EI]) as well as 

more exclusive event shapes like jet broadening [72Vl76j . 
Other event shape moments were considered at 0(q;^) in 
Ref. [IS] in the context of the dispersive model for the 
1/Q power corrections. 
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Appendix A: Theory parameter scan 

In this Appendix we describe the method we use to 
estimate uncertainties in our analysis. We will briefly re- 



view the profile functions and the theoretical parameters 
which determine the theory uncertainty. We will also de- 
scribe the scan over those parameters and the effects they 
have on the fit results. 

The profile functions used in Ref. [H], to which we re- 
fer for a more extensive description, are r-dependent fac- 
torization scales which allow us to smoothly interpolate 
between the theoretical constraints the hard, jet and soft 
scale must obey in different regions of the thrust distri- 
bution: 



1) peak: 

2) tail: 

3) far-tail: 



fJ'H ^ Q, 
^J■H = fJ-j 



\J Aqcd ' 



(Al) 



The factorization theorem derived for thrust in Ref. [S] 
is formally invariant under 0(1) changes of the profile 
function scales. The residual dependence on the choice 
of profile functions constitutes one part of the theoretical 
uncertainties and provides a method to estimate higher 
order perturbative corrections. We adopt a set of six 
parameters that can be varied in our theory error analysis 
which encode this residual freedom while still satisfying 
the constraints in Eq. (Al). 



For the profile function at the hard scale, we adopt 

m=eHQ, (A2) 

where eu is a free parameter which we vary from 1/2 to 
2 in our theory error analysis. 

For the soft profile function we use the form 




r)^ 



< T < ^1, 

ti<T< h, (A3) 
t2<r<\. 



Here, ti and t2 represent the borders between the peak, 
tail and far-tail regions, /ig is the value of /ig at r = 0. 
Since the thrust value where the peak region ends and the 
tail region begins is Q dependent, ti ~ l/Q, we define 
the Q-independent parameter ni by ti — ni/ {Q /I GeV). 
To ensure that hs{t) is a smooth function, the quadratic 
and linear forms are joined by demanding continuity of 
the function and its first derivative at t = ti and t = t2, 
which fixes 6 = 2 (/i/f — /io) 7(^2 — ii + 5) and d = [/Lto(i2 + 
\) — IJ-H ii] 7(^2 ~ ^1 + 5) ■ In our theory error analysis we 
vary the free parameters ni, t2 and /io. 

The profile function for the jet scale is determined by 
the natural relation between the hard, jet, and soft scales 

A^j(t) = (^l + ej(^-T)'^ y^^M^- (A4) 

The term involving the free 0(l)-parameter ej imple- 
ments a modification to this relation and vanishes in the 
multijet region where r — 1/2. We use a variation of ej 
to include the effect of such modifications in our estima- 
tion of the theoretical uncertainties. 
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FIG. 16: Impact on parameters of the M\ fit from variations of the best-fit values for asijnz) and f2i values in the ranges 
given in Table |VIIl] The dark shaded blue regions represent values of the parameters larger than their default values, the light 
shaded green regions where the parameters are smaller than their default values. 



In our theory error analysis we vary to account for 
our ignorance on the resummation of logarithms of r in 
the nonsingular corrections. We consider three possibili- 
ties 



Mns(T) { Hj{t), Us 0, 



(A5) 



The complete set of theoretical parameters and the 
their ranges of variation are summarized in Table |VIII[ 

Besides the parameters associated with the profile 
functions, the other theory parameters are T™^^, J3, 
S3, and €2.3. The cusp anomalous dimension at ©(a^), 
Y™^^ is estimated via Pade approximants and we assign a 
200% uncertainty to this approximation, js and S3 repre- 
sent the nonlogarithmic 3-loop term in the position-space 
hemisphere jet and soft functions, respectively. These 
two parameters and their variations are estimated via 
Pade approximations. The last two parameters 62 and £3 
allow us to include the statistical errors in the numeri- 
cal determination of the nonsingular distribution at two 
(from EVENT2 [771 EH]) and three (from EERAD3 [2j) 
loops, respectively. 

At each order we randomly scan the parameter space 
summarized in Table |VIII| with a uniform measure, ex- 
tracting 500 points. Each of the points in Fig. [6] is the 
result of the fit performed with a single choice of a point 
in the parameter space. The contour of the area in the 
as-2 rii plane covered by the fit results at each given or- 
der is fitted to an ellipse, which is interpreted as a 1-cr 



theoretical uncertainty. The ellipse is determined as fol- 
lows: in a first step we determine the outermost points 
on the ag-'i^li plane (defined by the outermost convex 
polygon). We then perform a fit to these points using a 
which is the square of the formula for an ellipse: 



Xcllips 



(A6) 



Here the sum is over the outermost points. The coordi- 
nates for the center of the ellipse, ag and fip, are fixed 
ahead of time to the average of the maximum and min- 
imum values of as{mz) and fii in the scan. We then 
minimize XciUpsc determine the parameters a, 6, c of 
the ellipse. 

One could further express the coefficients a and b by 



1 + Vl + 4c2 Aa2 Af72 
2Aa2 

1 + \/l + 4 c2 Aa2 
8Af72 



(A7) 



where Aa and Ail are just the half of the difference of 
the maximum and minimum values of as{mz) and fii, 
respectively, on the ellipse. Setting Aa and AJ7 to the 
corresponding values obtained from the fit points of the 
scan (i.e. the perturbative errors) the coefficients a and b 
can be fixed and only c remains as a free parameter. The 
minimization of Xoiiipso ^'i- (^6) gives almost iden- 



tical results regardless of whether or not Eqs. (A7| are 
imposed. 
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In Fig. 16 we vary a single parameter of Table VIII 
keeping all the others fixed at their respective default 
values, and we plot the change of as{mz) and fii as com- 
pared to the values obtained from the first moment thrust 
fit with the default setup. In the figure, the dark shaded 
blue area represents a variation where the parameter is 
larger than the default value, and the light shaded green 
one where the parameter is smaller. The largest uncer- 
tainty is associated with the variation of the hard scale, 
en- The value of as{mz) is similarly affected by the un- 



certainty of the profile function parameters, the statisti- 
cal error from the numerical determination of the 3-loop 
nonsingular distribution from EERAD3 [5], and by the 
parameter J3. It is rather insensitive to the variation of 
the 4-loop cusp anomalous dimension and the statistical 
error from the determination of the 2-loop nonsingular 
contribution to the thrust distribution. The value of fii 
is mainly sensitive to the profile function parameters and 
£3, but is quite insensitive to J3. 
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